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P2003-012526 filed on January 21, 2003; the entire contents of 
which are incorporated by reference herein. 

BACKGROUND OF THE INVENTION 
1. Field of the Invention 

The present invention relates to a computer 
implemented method for development profile simulation, a 
computer program product for controlling a computer system 
so as to simulate development profile, and a computer 
implemented method for mask pattern data correction and in 
particular to a simulation method to predict a resist pattern 
shape after exposure and development and a computer program 
for the simulation method and a mask pattern data correction 
method using the predicted pattern shape. 
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2. Description of the Related Art 

Generally, a series of processes for manufacturing 

semiconductor devices includes a lithography process to make 

a fine pattern on a substrate. In the lithography process, 

5 various parameters such as exposure conditions for a 

projection system and resist process conditions are present to 

provide a desired resist pattern. For example, the exposure 

conditions for the projection system include an aperture, a 

coherence factor, the shape of a light source and a pupil filter. 

10 The resist process conditions include the types and 
thicknesses of resist and anti-reflection films, a post-exposure 
baking (PEB) temperature, a PEB time, and development time. 
Furthermore, the parameters for the lithography process also 
include reticle conditions such as mask pattern size and mask 

15 type to provide the desired resist pattern. 

As semiconductor devices have continued to shrink in 
size over time, it is impossible to obtain the desired resist 
pattern with a desired depth of focus without optimizing the 
parameters described above, because of a narrow margin of the 

20 depth of focus. However, optimizing the parameters for the 
lithography process by actually manufacturing masks over and 
over requires a great deal of labor, time and costs. Therefore, 
in order to reduce the labor, time and costs, a computer 
simulation is employed to estimate the shape of a developed 

25 photosensitive resist. Various resist profile simulation models 
such as a string model, a ray tracing model and a cell model 
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having been proposed. For the existing resist profile 
simulation models, it is assumed that the dissolution rate of 
the photosensitive resist with the developer depend only on the 
exposure dose, the photosensitive character of the resist and 
5 the process condition, or only on the resist pattern shape at a 
dissolved portion in addition to the above described items. 

However, the dissolution rate of the photosensitive 
resist during the actual development is affected by the change 
in alkaline concentration. For example, since the resist is 

10 developed by the neutralization reaction between an acid unit 
in the resist and the alkaline in the developer, the alkaline 
concentration of the developer varies as the neutralization 
reaction progresses. Therefore, it is difficult to assume that 
the existing resist profile simulation models accurately 

15 predict the resist pattern shape and the resist pattern size. 

In Japanese Patent LaidOpen Publication No. Hei 
10*64792, a resist profile simulation model that takes into 
account the change in the alkaline concentration is proposed. 
According to the proposed resist profile simulation model, the 

20 dissolution rate is adjusted in accordance with either the 
spatial average value of a latent image intensity or a 
photosensitizer concentration on the assumption that the 
change in the alkaline concentration is proportional to the 
amount of dissolved resist around a target resist pattern. The 

25 proposed resist profile simulation model assumes that the 
changing ratio of the dissolution rate is proportional to the 
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amount of the dissolved resist around the target resist pattern, 
i.e., the change in the alkaline concentration. 

However, when the dissolution rate is calculated by an 
actual measurement, the changing ratio of the dissolution rate 
5 is not always proportional to the change in the alkaline 
concentration. Therefore, the existing resist profile simulation 
models may fail to converge on a calculation to optimize the 
parameters. Even when the parameters are optimized, the 
predicted resist pattern size obtained by the simulation may 

10 not be equal to the actual resist pattern size obtained by 
actual measurement. 

Further, when a developer is supplied to the resist, the 
developer flows in a direction in accordance with the flow 
velocity. In the case of a scan development method where a 

15 developer nozzle is scanned a wafer from one end to the other, 
the flow of the developer may occur during the development, 
after a puddle is formed. Since the alkaline concentration of 
the developer is reduced by the dissolution of the resist, the 
concentration of the developer varies on the concentration of 

20 the resist pattern. Therefore, the actual resist pattern size 
may differ from the predicted resist pattern size. Consequently, 
the existing resist profile simulation models also fail to 
predict the resist pattern size variation due to the flow of the 
developer. 

25 

SUMMARY OF THE INVENTION 
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An aspect of present invention inheres in a computer 
implemented method for development profile simulation 
according to an embodiment of the present invention including 
calculating optical intensities in a photosensitive resist, 
5 calculating a spatial average value of the optical intensities, 
reading a measured changing ratio of a dissolution rate of the 
photosensitive resist relating to an alkaline concentration 
changed by at least one of exposure dose on the photosensitive 
resist, a position in the thickness direction of the 
10 photosensitive resist and an alkaline concentration of 
developer for the photosensitive resist, obtaining a calculated 
dissolution rate by using the spatial average value and the 
measured changing ratio, and predicting a pattern shape of the 
photosensitive resist from the calculated dissolution rate. 
15 Another aspect of the present invention inheres in the 

computer implemented method for development profile 
simulation according to the embodiment of the present 
invention including calculating optical intensities in a 
photosensitive resist, calculating a spatial average value of 
20 the optical intensities, o btaining a changing ratio of a 
logarithm of a measured dissolution rate to an alkaline 
concentration of developer for the photosensitive resist or the 
changing ratio of the logarithm of the measured dissolution 
rate to a logarithm of the alkaline concentration of developer 
25 for the photosensitive resist, obtaining a calculated 
dissolution rate by using the spatial average value and the 
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calculated changing ratio of the logarithm of the measured 
dissolution rate to alkaline concentration of a developer or the 
calculated changing ratio of the logarithm of the measured 
dissolution rate to the logarithm of an alkaline concentration 
of the developer, and predicting a pattern shape of the 
photosensitive resist by using the calculated dissolution rate. 

Yet another aspect of the present invention inheres in 
the computer implemented method for development profile 
simulation according to the embodiment of the present 
invention including calculating optical intensities in a 
photosensitive resist, calculating a spatial average value of 
the optical intensities in an area where the photosensitive 
resist dissolves in initial stage of development, reading a 
measured changing ratio of a dissolution rate of the 
photosensitive resist relating to an alkaline concentration, 
obtaining a calculated dissolution rate by using the spatial 
average value and the measured changing ratio, and predicting 
a pattern shape of the photosensitive resist by using the 
calculated dissolution rate. 

Yet another aspect of the present invention inheres in 
the computer implemented method for development profile 
simulation according to the embodiment of the present 
invention including calculating optical intensities in a target 
position to predict a pattern shape of a photosensitive resist 
and reference positions moving along development time, 
calculating spatial average values of the optical intensities in 



the reference positions, reading a measured changing ratio of a 
dissolution rate of the photosensitive resist relating to the 
alkaline concentration, obtaining calculated dissolution rates 
by using the spatial average values in the reference positions 
5 and the measured changing ratio, and predicting the pattern 
shape of the photosensitive resist in the target position by 
using the calculated dissolution rates and the optical 
intensities in the target position. 

Yet another aspect of the present invention inheres in 

10 a computer program product for controlling a computer 
system so as to simulate development profile, the computer 
program product according to the embodiment of the present 
invention including instructions configured to calculate 
optical intensities in a photosensitive resist within the 

15 computer system, instructions configured to calculate a 
spatial average value of the optical intensities within the 
computer system, instructions configured to read a measured 
changing ratio of a dissolution rate of the photosensitive resist 
relating to an alkaline concentration changed by at least one of 

20 exposure dose on the photosensitive resist, a position in the 
thickness direction of the photosensitive resist and an alkaline 
concentration of developer for the photosensitive resist within 
the computer system, instructions configured to obtain a 
calculated dissolution rate by using the spatial average value 

25 and the measured changing ratio within the computer system, 
and instructions configured to predict a pattern shape of the 
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photosensitive resist from the calculated dissolution rate 
within the computer system. 

Yet another aspect of the present invention inheres in the 
computer program product for controlling a computer system 
so as to simulate development profile, the computer program 
product according to the embodiment of the present invention 
including instructions configured to calculate optical 
intensities in a photosensitive resist within the computer 
system, instruction configured to calculate a spatial average 
value of the optical intensities within the computer system, 
instruction configured to obtain a changing ratio of a 
logarithm of a measured dissolution rate to an alkaline 
concentration of developer for the photosensitive resist or the 
changing ratio of the logarithm of the measured dissolution 
rate to a logarithm of the alkaline concentration of the 
developer for the photosensitive resist within the computer 
system, instruction configured to obtain a calculated 
dissolution rate by using the spatial average value and the 
calculated changing ratio of the logarithm of the measured 
dissolution rate to alkaline concentration of a developer or the 
calculated changing ratio of the logarithm of the measured 
dissolution rate to the logarithm of the alkaline concentration 
of the developer within the computer system, and instruction 
configured to predict a pattern shape of the photosensitive 
resist by using the calculated dissolution rate within the 
computer system. 
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Yet another aspect of the present invention inheres in the 
computer program product for controlling a computer system 
so as to simulate development profile, the computer program 
product according to the embodiment of the present invention 
5 including instruction configured to calculate optical 
intensities in a photosensitive resist within the computer 
system, instruction configured to calculate a spatial average 
value of optical intensities in an area where the photosensitive 
resist dissolves in initial stage of development within the 

10 computer system, instruction configured to read a measured 
changing ratio of a dissolution rate of the photosensitive resist 
relating to an alkaline concentration within the computer 
system, instruction configured to obtain a calculated 
dissolution rate by using the spatial average value and the 

15 measured changing ratio within the computer system, and 
instruction configured to predict a pattern shape of the 
photosensitive resist by using the calculated dissolution rate 
within the computer system. 

Yet another aspect of the present invention inheres in the 

20 computer program product for controlling a computer system 
so as to simulate development profile, the computer program 
product according to the embodiment of the present invention 
including instruction configured to calculate optical 
intensities in a target position to predict a pattern shape of a 

25 photosensitive resist and reference positions moving along 
development time within the computer system, instruction 
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configured to calculate spatial average values of the optical 
intensities in the reference positions within the computer 
system, instruction configured to read a measured changing 
ratio of a dissolution rate of the photosensitive resist relating 
5 to the alkaline concentration within the computer system, 
instruction configured to obtain calculated dissolution rates 
by using the spatial average values in the reference positions 
and the measured changing ratio within the computer system, 
and instruction configured to predict the pattern shape of the 

10 photosensitive resist in the target position by using the 
calculated dissolution rates and the optical intensities in the 
target position within the computer system. 

Yet another aspect of the present invention inheres in a 
computer implemented method for mask pattern data 

15 correction according to the embodiment of the present 
invention including reading a designed pattern data in a 
photosensitive resist, a mask pattern data, and a measured 
changing ratio of a dissolution rate of the photosensitive resist 
relating to an alkaline concentration changed by at least one of 

20 exposure dose on the photosensitive resist, a position in the 
thickness direction of the photosensitive resist and an alkaline 
concentration of developer for the photosensitive resist, 
calculating optical intensities in the photosensitive resist by 
using the mask pattern data, calculating a spatial average 

25 value of the optical intensities, obtaining a calculated 
dissolution rate by using the spatial average value and the 
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measured changing ratio, predicting a pattern shape of the 
photosensitive resist from the calculated dissolution rate, and 
optimizing the mask pattern data so as to make the 
calculated pattern shape similar to the designed pattern data 
5 in the photosensitive resist. 

Yet another aspect of the present invention inheres in the 
computer implemented method for mask pattern data 
correction according to the embodiment of the present 
invention including obtaining a designed pattern data in a 

10 photosensitive resist, a mask pattern data, and a calculated 
changing ratio of a calculated changing ratio of a logarithm of 
a measured dissolution rate to an alkaline concentration of 
developer for the photosensitive resist or the calculated 
changing ratio of the logarithm of the measured dissolution 

15 rate to a logarithm of the alkaline concentration of developer 
for the photosensitive resist, calculating optical intensities in 
the photosensitive resist by using the mask pattern data, 
calculating a spatial average value of the optical intensities, 
obtaining a calculated dissolution rate by using the spatial 

20 average value and the measured changing ratio, predicting a 
pattern shape of the photosensitive resist from the calculated 
dissolution rate, and optimizing the mask pattern data so as to 
make the calculated pattern shape similar to the designed 
pattern data in the photosensitive resist. 

25 Yet another aspect of the present invention inheres in the 
computer implemented method for mask pattern data 



correction according to the embodiment of the present 
invention including reading a designed pattern data in a 
photosensitive resist, a mask pattern data, and a measured 
changing ratio of a dissolution rate of the photosensitive resist 
5 relating to an alkaline concentration, calculating optical 
intensities in the photosensitive resist by using the mask 
pattern data, calculating a spatial average value of optical 
intensities in an area where the photosensitive resist 
dissolves in initial stage of development, obtaining a 

10 calculated dissolution rate by using the spatial average value 
and the measured changing ratio, predicting a pattern shape of 
the photosensitive resist from the calculated dissolution rate, 
and optimizing the mask pattern data so as to make the 
calculated pattern shape similar to the designed pattern data 

15 in the photosensitive resist. 

Yet another aspect of the present invention inheres in the 
computer implemented method for mask pattern data 
correction according to the embodiment of the present 
invention including reading a designed pattern data in a 

20 photosensitive resist, a mask pattern data, and a measured 
changing ratio of a dissolution rate of the photosensitive resist 
relating to the spatial average value, calculating optical 
intensities in a target position of the photosensitive resist and 
reference positions moving along development time by using 

25 the mask pattern data, calculating a spatial average value of 
the optical intensities in the reference positions, obtaining a 
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calculated dissolution rate by using the spatial average value 
in the reference positions and the measured changing ratio, 
predicting a pattern shape of the photosensitive resist by 
using the calculated dissolution rate and the optical 
intensities in the target position, and optimizing the mask 
pattern data so as to make the calculated pattern shape 
similar to the designed pattern data in the photosensitive 
resist. 

BRIEF DESCRIPTION OF DRAWINGS 

Fig.l is a graph that explains a way to obtain a 
measured changing ratio of a dissolution rate "g"> 

Fig. 2 is a graph that illustrates the relation between a 
changing ratio of the dissolution rate and a resist position in 
the thickness direction; 

Fig.3A is a graph that illustrates the dependence of the 
changing ratio of the dissolution rate on exposure dose in a top 
position of the resist and Fig.3B is a graph that illustrates the 
dependence of the changing ratio of the dissolution rate on 
exposure dose in a middle position of the resist; 

Fig. 4 is a diagram of an apparatus for a profile 
simulation in accordance with a first embodiment of the 
present invention; 

Fig. 5 is a flowchart of a method for optimizing 
development parameters in accordance with the first 
embodiment of the present invention; 
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Fig. 6 is a flowchart of a method for development profile 
simulation in accordance with the first embodiment of the 
present invention; 

Fig. 7 is a plan view of a layout pattern in accordance 
with the first embodiment of the present invention; 

Fig.8A is a graph that illustrates a actual hole diameter 
relating to position and Fig.8B is a graph that illustrates a 
calculated hole diameter relating to the position; 

Fig. 9 is a block diagram of the apparatus for the profile 
simulation; 

Fig. 10 is a flowchart of a method for optimizing 
development parameters and photoreaction parameters in 
accordance with the first embodiment of the present invention; 

Fig. 11 is a second flowchart of a method for optimizing 
development parameters and photoreaction parameters in 
accordance with the first embodiment of the present invention; 

Fig. 12 is a flowchart of the method for development 
profile simulation in accordance with the first embodiment of 
the present invention; 

Fig. 13 is a sectional diagram that illustrates a flow of 
the developer on the resist in accordance with the first 
embodiment of the present invention. 

Fig. 14 is a plan view of a layout pattern in accordance 
with a second embodiment of the present invention; 

Fig.l5A is a graph that illustrates a hole diameter with 
non-corrected photo mask relating to a position and Fig.lSB is 
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a graph that illustrates the hole diameter with a corrected 
photo mask relating to the position; 

Fig. 16 is a block diagram of an apparatus for mask 
pattern data correction in accordance with a third embodiment 
5 of the present invention; and 

Fig. 17 is a flowchart of a method for mask pattern data 
correction in accordance with the third embodiment of the 
present invention; 

10 DETAILED DESCRIPTION OF THE INVENTION 

Various embodiments of the present invention will be 
described with reference to the accompanying drawings. It is 
to be noted that the same or similar reference numerals are 
applied to the same or similar parts and elements throughout 

15 the drawings, and the description of the same or similar parts 
and elements will be omitted or simplified. 

(Development Model) 

Before describing a preferred embodiment of the present 
20 invention, a development model for a resist development 
process is described below. The development of the resist is a 
neutralization reaction between an acid unit in the resist and 
the alkaline in a developer. In the development process, the 
neutralization reaction between the resist and the developer 
25 progresses, and the alkaline concentration of the developer 
decreases as the resist dissolves. The dissolution rate for a 
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highly exposed area of a positive resist or for a low exposed 
area of a negative resist is very high, for example, 1,000 
nm/sec to 10,000 nm/sec. Therefore, such area of the resist 
dissolves at the initial development stage. 
5 After the initial development stage, a low exposure area 

such as a pattern edge portion dissolves. Since the resist has 
already dissolved partially at the initial stage, it is assumed 
that the developer of a low alkaline concentration makes the 
low exposure area of the pattern edge dissolve and 

10 consequently the final resist pattern size is determined. 

In many cases, the relation between the alkaline 
concentration of the developer and the dissolution rate is not 
proportional to each other, and as the alkaline concentration 
increases, the dissolution rate also increases exponentially. 

15 Therefore, it is appropriate that the relation between the 
alkaline concentration and the dissolution rate is given by the 
approximate equation (l) or (2)' 

d{log R(c)} /5c = const. (l) 

20 

5{log R(c)} Id (log c) = const. (2) 

Here, c is alkaline concentration of the developer, and 
R(c) is the dissolution rate of the resist. The changing ratio of 
25 the logarithm of the dissolution rate R(c) of a photosensitive 
resist to the alkaline concentration of the developer as shown 
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in equation (3) is used as an example of a changing ratio "g" 
that shows the changing ratio of the resist dissolution rate 
relating to the alkaline concentration of the developer. When 
the initial alkaline concentration of the supplied developer is 
co and the initial dissolution rate is R(co), a measured 
dissolution rate R(c) is given by the equation (4): 

a{log R(c)} Idc = g (3) 

log R(c) = g(c-co) + log R(c 0 ) (4) 

The quantity of the reacting acid unit of the resist is 
proportional to the amount of the dissolved resist around the 
resist pattern. The change of the alkaline concentration is 
obtained by subtracting the amount of alkaline used by the 
neutralization reaction from the initial alkaline concentration 
co. At a predetermined position (x,y) in the resist, the 
amount of dissolved resist around the resist pattern, i.e., the 
amount of alkaline used by the neutralization reaction, is 
given by the equation (5). The spatial average value G(x, y) is 
given by the equation (6). 

c = co ' a/G(x, y) (5) 
G(x,y) 

= |Jdx ) dy , I(x,y)exp[{-(x-x , ) 2 + (y-y , ) 2 }/k2] ( 6 ) 
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Here, a* is a constant. G(x,y) is the spatial average value 
of an aerial image intensity and I(x,y) is an aerial image 
intensity I(x,y). 

The relation between the dissolution rate R(x,y) and the 
spatial average value G(x,y) in the predetermined position 
(x,y) is given by the equation (7) based on the equations (4) and 
(5): 

log R(x, y) = gocG(x, y) + log R(c 0 ) (7) 

Here, a is -<x\ and get is the changing ratio of the 
dissolution rate to the spatial average value of the aerial 
image intensity. 

In this case, the aerial image intensity I(x,y) is 
employed to obtain G(x,y), however, I(x,y) may be replaced 
with M(x,y) concentration of photore action products generated 
during exposure and PEB. 

The changing ratio of the dissolution rate "g" shown in 
the equation (3) is obtained through measured measurement. 
For example, by changing the alkaline concentration of the 
supplied developer, it is possible to measure the dissolution 
rate corresponding to the alkaline concentration with a 
dissolution rate monitor (DRM). In Fig. 1, the relation between 
the alkaline concentration "c" and the logarithm of the 
dissolution rate R(x,y) is shown and the results obtained by 
the DRM are plotted. By drawing an approximation curve 22 
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through a plurality of results 21a to 21c, it is possible to 
obtain the measured changing ratio "g" (Gradient) of the 
dissolution rate. 

The behavior of the changing ratio of the dissolution 
rate "g" measured by the DRM depends on the resist type. For 
example, there is a resist in which the changing ratio of the 
dissolution rate "g" hardly depends on the exposure dose and 
resist position in the thickness direction. In this case, the 
changing ratio "g" obtained by plotting the dissolution rates as 
shown in Fig.l shows almost the same value although exposure 
dose and resist positions in the thickness direction is changed. 
The approximation curve 22 in Fig. 1 is linear through out the 
overall alkaline concentration. Therefore, it is possible to 
modulate the dissolution rate R(x,y) under the assumption 
that the changing ratio of the dissolution rate "g" is constant. 

On the other hand, there is a different resist (resist A) 
in which the changing ratio of the dissolution rate "g" depends 
on the resist position in the thickness direction, as shown in 
Fig. 2. In Fig. 2, the changing ratio of the dissolution rate "g" 
is the average value in terms of exposure dose. Along the 
horizontal axis, 0 nm indicates the bottom of the resist, and 
400 nm indicates the top of the resist. As shown in Fig. 2, the 
changing ratio of the dissolution rate "g" is constant in the 
interior of the resist, while the changing ratio of the 
dissolution rate "g" in the upper and lower portions of the 
resist varies. Therefore, the changing ratio of the dissolution 
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rate "g" is defined for the upper, the interior and the lower 
portions of the resist, respectively. In a position "p" inside the 
resist through the range of 32 nm to 205 nm of the thickness 
direction, g=23.6 which is an average value defined as the 
5 changing ratio of the dissolution rate. In the upper and lower 
portions of the resist, changing ratio of the dissolution rate "g" 
is given by the linear approximate equations (8). and (9). In 
the method for development profile simulation described below, 
a table is formed containing the changing ratio of the 
10 dissolution rates "g" (Gradient) for the individual resist, 
positions in the thickness direction, and the changing ratio of 
the dissolution rate "g" varies in accordance with the resist 
position in the thickness direction. 

15 g = 0.1673 p + 18.08 (p = 0~32) (8) 

g = 0.023 p + 18.61 (p = 205~ 379) (9) 

Furthermore, there is another resist in which the 
20 changing ratio of the dissolution rate "g" slightly depends on 
the exposure dose or the alkaline concentration of the 
developer. In the case where the changing ratio of the 
dissolution rate "g" depends on the exposure dose, it is 
possible to define the changing ratio of the dissolution rate "g" 
25 by replacing the horizontal axis in Fig. 2 with the exposure 
dose. Therefore, the changing ratio of the dissolution rate "g" 
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is defined respectively in areas where the changing ratio "g" 
does not depend on the exposure does and in areas where the 
changing ratio "g" depends on the exposure dose. 

In the case where the changing ratio of the dissolution 
5 rate "g" depends on the alkaline concentration "c" of the 
developer, the inclination slope of the approximation curve 22 
shown in Fig. 1 (the changing ratio of the dissolution rate "g") 
is not constant. In this case, a range of the alkaline 
concentration "c" is divided into a plurality of segments such 

10 as low, middle and high. The inclination slope of the 
approximation curve 22 for each segment is defined as the 
changing ratio of the dissolution rate "g". Thereafter, a table is 
formed containing the changing ratio of the dissolution rate 
"g" relating to the alkaline concentration "c". It is assumed 

15 that the amount of change of the alkaline concentration c due 
to the development process is proportional to the spatial 
average value G(x, y) of the aerial image intensity or the 
concentration of photoreaction products. Therefore, in the 
method for development profile simulation described below, 

20 the dependency on the alkaline concentration of the changing 
ratio of the dissolution rate "g" may be defined relative to the 
spatial average value G(x, y). In this case, a range of the 
spatial average value G(x, y) may be divided into a plurality of 
segments so that each segment corresponds to the alkaline 

25 concentration. For each segment, the table containing the 
changing ratio of the dissolution rate "g" relating to the 
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spatial average value G(x, y) may be formed. 

In addition, there is another resist (resist B) in which 
the changing ratio of the dissolution rate "g" is drastically 
changed relative to the alkaline concentration of the developer, 
the exposure dose and the resist position in the thickness 
direction, as shown in Figs. 3A and 3B. 

In Figs. 3A and 3B, the changing ratio of the dissolution 
rate "g" is represented relative to the exposure dose. The 
alkaline concentration is divided, i.e., 0.21N to 0.24N and 
0.24N to 0.27N. It is apparent from Figs. 3A and 3B that the 
changing ratio of the dissolution rate "g" greatly depends on 
the resist position in the thickness direction, the exposure 
dose and the alkaline concentration of the developer. Based 
on the results shown in Figs. 3A and 3B, a table is formed 
containing the ratio of the logarithm of the dissolution rate to 
the alkaline concentration as the measured value, as the 
function of the resist position in the thickness direction, the 
exposure dose, and the alkaline concentration. Alternatively, 
the measured values of the changing ratio of the dissolution 
rate "g"s shown in Figs. 3A and 3B may be entered in the table. 
However, considering the measurement error included in the 
measured value, an average values should be entered as the 
constant value where a linear approximate is obtained. For 
example, in Fig. 3A, it is possible to use a constant value for 
the range where the exposure dose is less than 5.7 mJ and the 
exposure dose is more than 6.0 mJ while the alkaline 
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concentration is from 0.24N to 0.27N. Further, in Fig. 3A, the 
linear approximation is conducted for the range where the 
exposure dose is 5.7 mJ to 6.0 mj while the alkaline 
concentration is 0.24N to 0.27N. 

The above explanation describes a method to obtain the 
changing ratio of the dissolution rate "g" calculated by the 
equation (3). However, the function approximating the relation 
between the alkaline concentration c and the dissolution rate 
R(x, y) is not limited to the equation (3). Even though the 
relation between the alkaline concentration c and the 
dissolution rate R(x, y) may be represented by an equation 
other than the equation (3), it is also expected 
phenomenologically that the changing ratio of the dissolution 
rate may change in accordance with the resist position in the 
thickness direction, the exposure dose and the alkaline 
concentration. 

Segregation of a photo-acid generator (PAG) and 
evaporation of the acid may change the dependence of the 
changing ratio of the dissolution rate "g" on the resist position 
in the thickness direction. Some resist is designed so that the 
PAG is segregated on the surface of a resist in order to form a 
rectangular resist pattern. For such kind of the resist, since a 
dissolution inhibitation effect due to the PAG differs between 
the surface and the interior of the resist, the changing ratio of 
the dissolution rate "g" may depend on the resist position in 
the thickness direction. It is also well known that, since the 
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acid is evaporated from the surface of the resist, the reduction 
of acid on the surface of the resist decreases the dissolution 
rate of the surface of the resist. Therefore, the changing 
ratio of the dissolution rate "g" may depend on the resist 
5 position in the thickness direction. 

The variation in the changing ratio of the dissolution 
rate "g" may also be caused by the fact the resist contains a 
mixture of different types of resins. In order to maintain the 
shape of the resist and to reduce roughness, a resist containing 

10 a different molecular weight of resin or a different type of 
resin is generally used. The dissolution probability of a resin 
is varied in accordance with the number of protecting groups 
decomposed by exposure. Therefore, in the case where 
different resins are included in the resist, the dissolution rate 

15 of the resin may differ on the exposure dose. Hence, in the 
case where the changing ratio of the dissolution rate "g" 
differs for the individual resins, the changing ratio of the 
dissolution rate "g" of the resist also changes based on the 
exposure dose. Further, the dissolution probability of the 

20 resin in the resist also may vary in accordance with the initial 
alkaline concentration Co. Therefore, in the case where 
different resins are present in the resist, the changing ratio of 
the dissolution rate "g" may change depending on the initial 
alkaline concentration Co. 

25 Since the dissolution behavior of the resist has not been 

completely investigated, it is impossible at present time to 
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explain all phenomena where the changing ratio of the 
dissolution rate "g" corresponding to the change of the alkaline 
concentration c may change. However, even when the 
function approximating the relation between the alkaline 
concentration c and the dissolution rate R(x, y) is changed 
from equation (3) to another equation (e.g., equation (2)), it is 
appropriate to change the changing ratio of the dissolution 
rate u g" corresponding to the change of the alkaline 
concentration as the function of resist position in the 
thickness direction, the exposure dose, and the alkaline 
concentration. 

By employing the above described changing ratio of the 
dissolution rate "g" corresponding to the alkaline 
concentration c, a method for simulating the resist pattern 
shape is described below. 

(First Embodiment) 

As shown in Fig. 4, an apparatus for development profile 
simulation according to a first embodiment of the present 
invention includes a processing unit 1 having a function for 
performing the method for development profile simulation, a 
data memory 2 and a program memory 3 connected to the 
processing unit 1. The data memory 2 stores various data, 
such as an exposure condition, mask pattern data and a 
development parameter 12 which are used to operate the 
processing unit 1. The program memory 3 stores various 
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types of application software, middleware and an operating 
system that are used to operate the processing unit 1. 

The processing unit 1 includes optical intensities 
calculator configured to calculate optical intensities of a 
5 photosensitive resist, a spatial average value calculator 8 
configured to calculate a spatial average value of the optical 
intensities, a changing ratio of the dissolution rate reader 9 
configured to read an exposure condition, mask pattern data, 
and a measured changing ratio of the dissolution rate of the 

10 resist relating to an alkaline concentration, a dissolution rate 
calculator 10 configured to obtain a calculated dissolution rate 
using the spatial average value and the measured changing 
ratio of the dissolution rate, and a pattern shape calculator 11 
configured to predict the pattern shape of the resist by using 

15 the calculated dissolution rate. 

The "optical intensities" includes an aerial image 
intensity and concentration of photoreaction products. In the 
first embodiment, the aerial image intensity is employed as 
the optical intensities. Therefore, the "optical intensities 

20 calculator" is an aerial image intensity calculator 7 and the 
spatial average value calculator 8 calculates the spatial 
average value of the aerial image intensity. Further, in the 
first embodiment, the ratio of the logarithm of the measured 
dissolution rate of the resist to the alkaline concentration is 

25 employed as the changing ratio of the dissolution rate "g" of 
the resist to the alkaline concentration. 
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The processing unit 1 constitutes one part of the central 
processing unit (CPU) of an ordinary computer system. The 
aerial image intensity calculator 7, the spatial average value 
calculator 8, the changing ratio of the dissolution rate reader 9, 
the dissolution rate calculator 10 and the pattern shape 
calculator 11 may be constituted by special hardware, or may 
be software providing a substantially equivalent function by 
using the CPU. Although not shown in figure, the apparatus 
for development profile simulation includes a computer 
readable program to perform the method for development 
profile simulation and a main storage device, such as a random 
access memory (RAM) to temporarily store data processed by 
the processing unit 1, for example. The data memory 2 and the 
program memory 3 correspond to auxiliary storage unit such as 
read-only memories (ROMs) or magnetic disk drives. The 
processing unit 1 is connected to an input unit 6 for receiving 
data and instructions entered by an operator through an I/O 
controller 4, and an output unit 5 for outputting the pattern 
shape of the resist. A keyboard, a mouse, a light pen and a 
flexible disk drive are available for the input unit 6, and a 
display device and a printer are available for the output unit 
5. 

With reference next to Figs. 5 and 6, the method for 
development profile simulation according to the first 
embodiment of the present invention is described below. In 
Fig. 5, a method for optimizing development parameters "k" 
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and "a" used for the profile simulation is represented. The 
development parameters "k" and "a" are used in the equations 
(6) and (7). 

(a) In step S101, N resist patterns are formed, and the N 
resist pattern sizes are measured. An measured pattern size 
Wi,o (i = 1, 2, 3, . . ., N) of the i-th resist pattern is stored in the 
data memory 2 shown in Fig. 4. The number of resist patterns 
(N) should be larger than the number of parameters to be 
optimized. In the first embodiment, N is set to three or more 
since two development parameters "a" and "k" are optimized. 

(b) In step S102, the measured initial dissolution rate 
R(co) of the developer used for pattern formation of the resist, 
and the measured dissolution rate R(c) during the development 
are measured. For example, a set of developers having a 
different alkaline concentration "c" are prepared, and the 
plural measured dissolution rate R(c) at the alkaline 
concentration "c" are measured. 

(c) In step S103, as described with reference to Figs. 1 to 
3, the ratio of the logarithm of the measured dissolution rate 
R(c) to the alkaline concentration c is calculated by using the 
measured dissolution rate R(c) as the measured value of the 
changing ratio of the dissolution rate "g". Thereafter, a table is 
formed containing the measured value of the changing ratio of 
the dissolution rate "g". The measured value of the changing 
ratio of the dissolution rate "g" is either completely a constant 
value, or differs in accordance with the alkaline concentration 
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of the developer, the resist position in the thickness direction 
and the exposure dose. Further, in case that it is possible to 
obtain the approximate function of the logarithm of the 
measured dissolution rate R(c) relating to the alkaline 
concentration "c" of the developer, the measured value of the 
changing ratio of the dissolution rate "g" is obtained from the 
fitting function. Since the measured initial dissolution rate 
R(co) depends on the exposure dose, the measured value R(co) 
is obtained as a table or the fitting function. The obtained 
measured value of the changing ratio of the dissolution rate 
"g" is stored in the data memory 2. 

(d) In step S104, the exposure condition used for the 
measurement, the condition for the mask pattern, such as the 
pattern arrangement and the mask size, a threshold value 6 
that is described below and used as a reference for optimizing 
the development parameter, and the number of parameter sets 
"M" to be employed are entered in the processing unit 1. The 
conditions may be entered by reading data from the data 
memory 2, or by an operator using the input unit 6. The 
"exposure condition" used for the actual measurement includes 
the condition such as NA and a and the exposure dose. The 
"parameter set" is a combination of the development 
parameters "a" and "k" to be optimized. The parameter set 
containing the first parameter set "on" and "ki" to the M-th 
parameter set "ocm" and "kivi" is prepared. 

(e) In step S105, the aerial image intensity calculator 7 
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calculates an aerial image intensity I(x, y) by using the 
conditions of exposure and the mask patterns used for the 
actual measurement. Thereafter, in step S106, the spatial 
average value calculator 8 calculates a spatial average value 
5 G(x, y) of the aerial image intensity by the equation (6) with 
the calculated aerial image intensity I(x, y). In step S107, 
the dissolution rate calculator 10 reads the measured changing 
ratio of the dissolution rate from the data memory 2. 

(f) In step S108, the dissolution rate calculator 10 reads 
10 the n-th parameter set "ot n " and "k n " from the data memory 2. 

For the first calculation, n is defined as 1. In step S109, the 
dissolution rate calculator 10 calculates a calculated 
dissolution rate R(x, y) by using the spatial average value G(x, 
y) and the measured value of the changing ratio of the 

15 dissolution rate "g". Thereafter, in step S110, the pattern 
shape calculator 11 estimates the resist pattern shape by using 
the calculated dissolution rate R(x, y) and the aerial image 
intensity I(x, y) obtained in step S105. In step Sill, the 
pattern shape calculator 11 obtains a calculated pattern size 

20 Wi, n (i = 1, 2, 3, . . ., N) in a target position of the i'th resist 
pattern that corresponding to the measured pattern size Wi,o 
stored in the data memory 2 at step S101. 

(g) In step S112, the pattern shape calculator 11 judges 
the similarity between the calculated pattern size Wi, n and the 

25 measured pattern size Wi,o based on the threshold value 8. 
Specifically, if the sum value of the squared difference between 
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the measured pattern size Wi,o and the calculated pattern size 
Wi, n for the N resist patterns represented as Z | W i>n - Wi, 0 I 2 is 
smaller than the squared threshold value 5 2 , step S113 is the 
next procedure. When Z I Wi, n - W it0 I 2 is equal to or larger than 
the squared threshold value 6 2 , step S114 is the next procedure. 
It should be noted that 8 = (Z|W i( i - W i>0 | 2 ) 1/2 is employed as a 
threshold value for the first calculation (n = l). 

(h) In step S113, the pattern shape calculator 11 defines 
5 = (E | Wi, n " Wi,ol 2 ) 1/2 as a new threshold. In this case, "a n " 
and "k n " are stored in data memory 2 as a temporary parameter. 
In step S114, the pattern shape calculator 11 determines 
whether "n" is equal to "M". In this case, the n-th parameter 
set "a n " and "k n " is employed. When "n" is equal to "M", the 
processing in Fig. 5 is terminated. When "n" is not equal to 
"M", it is assumed that a parameter set to be calculated 
remains and step S115 is the next procedure. 

(i) In step S115, the pattern shape calculator 11 
prepares the n+l-th parameter set "ct n +i" and "k n + i", and 
performs the processes from steps S108 to S114 using the new 
parameter set "oc n +i" and "k n +i". The processes from steps 
S108 to S114 are repeated until "n" is equal to "M". 

As described above, by searching for the parameter 
minimizing the Z | Wi, n ' Wj,o I 2 which is sum value of the 
squared difference between the measured pattern size Wi.o and 
the calculated pattern size Wj, n , the development parameters 
"a" and "k" are optimized. 
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With reference next to Fig. 6, the method for 
development profile simulation method using the optimized 
development parameters "a" and "k" according to the first 
embodiment is described below. 

(A) In step S121, the exposure condition for the resist 
pattern for which the size is to be predicted, the mask pattern 
conditions such as the pattern arrangement and the mask size, 
the measured initial dissolution rate R(co) , the measured 
changing ratio of the dissolution rate "g", and the development 
parameters "a" and "k" stored in the data memory 2 are 
entered in the aerial image intensity calculator 7. The 
conditions may be entered by an operator using the input unit 
6. The exposure condition includes the lighting condition such 
as NA and a and the exposure dose. The development 
parameters "a" and "k" are the optimized parameters obtained 
through the processing in Fig. 5. 

(B) In step S122, the aerial image intensity calculator 7 
calculates the aerial image intensity I(x, y) for a 
predetermined position (x, y) of the resist with a set condition 
such as exposure and the mask pattern. In step S123, the 
average spatial value calculator 8 calculates the average 
spatial value G(x, y) of the aerial image intensity by the 
equation (6) with the calculated aerial image intensity I(x, y). 

(C) In step S124, the dissolution rate calculator 10 
obtains the calculated dissolution rate R(x, y) for the 

^predetermined position of the resist by using the average 
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spatial value G(x, y), the measured initial dissolution rate 
R(co) , the measured changing ratio of the dissolution rate "g", 
and the optimized development parameters "a" and "k". In 
step S125, the pattern shape calculator 11 predicts the pattern 
shape of the resist by using the calculated dissolution rate R(x, 
y) and the aerial image intensity I(x, y) obtained in step S122. 
Finally, in step S126, the pattern shape calculator 11 predicts 
the calculated pattern size W for the target position. 

In the above described method for development profile 
simulation, the dissolution rate in the predetermined position 
of the resist varies in accordance with the average spatial 
value of the optical intensities at the predetermined position 
of the resist. The measured changing ratio of the dissolution 
rate used in the method for development profile simulation is 
defined in accordance with at least one of the exposure dose on 
the resist, the resist position in the thickness direction, and 
the alkaline concentration of the developer. Therefore, the 
changing ratio of the dissolution rate "g" ocrelating to the 
average spatial value of the aerial image intensity is also 
defined in accordance with at least one of the exposure dose on 
the resist, the resist position in the thickness direction, and 
the alkaline concentration of the developer. Consequently, 
the method for development profile simulation reflecting the 
actual behavior of the development is realized. Further, the 
method for development profile simulation according to the 
first embodiment makes it possible to accurately predict the 
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shape and the size of the resist pattern. 

In the method shown in Fig. 5, the threshold value 5 is 
updated in order to minimize I I Wi f „ - W i>0 I 2 as the result of 
the squared difference between the measured pattern size Wi,o 
5 and the calculated pattern size Wi, n for each pattern. However, 
so long that the £ | Wi, n - Wj.o I 2 is minimized, a method is not 
limited to that shown in Fig. 5. 

The above explanation describes the the method for 
development profile simulation using the aerial image 

10 intensity I(x, y) as the optical intensities. However, the first 
embodiment is not limited to the simulation using the aerial 
image intensity. It is also possible to use "the concentration 
M(x,y) of the photoreaction products" reacted during exposure 
and a PEB as the optical intensities instead of the aerial image 

15 intensity I(x, y). Therefore, the method for development profile 
simulation using the concentration M(x,y) of the photoreaction 
products is described below. 

With reference now to Fig. 9, an apparatus for 
development profile simulation according to the first 
20 embodiment of the present invention includes the processing 
unit 1, the data memory 2 and the program memory 3 
connected to the processing unit 1. The data memory 2 stores 
various data, such as the exposure condition, mask pattern 
data, the development parameter 12 and a photoreaction 
25 parameter 32 which are used to operate the processing unit 1. 

The processing unit 1 includes an photoreaction 
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products concentration calculator 31, a spatial average value 
calculator 8 configured to calculate a spatial average value of 
the concentration M(x,y) of the photoreaction products, the 
changing ratio of the dissolution rate reader 9 configured to 
read the ratio of the dissolution rate of the resist to the 
alkaline concentration of the developer which changes relating 
to at least one of the exposure dose, the resist position in the 
thickness position and the alkaline concentration of the 
developer, a dissolution rate calculator 10 configured to obtain 
a calculated dissolution rate using the spatial average value 
and the measured changing ratio of the dissolution rate, and a 
pattern shape calculator 11 configured to predict the pattern 
shape of the resist by using the calculated dissolution rate. 

The processing unit 1 constitutes one part of the central 
processing unit (CPU) of an ordinary computer system. The 
photoreaction product concentration calculator 31, the spatial 
average value calculator 8, the changing ratio of the 
dissolution rate reader 9, the dissolution rate calculator 10 
and the pattern shape calculator 11 may be constituted by 
special hardware, or may be software providing a substantially 
equivalent function by using the CPU. The processing unit 1 is 
connected to the input unit 6 for receiving data and 
instructions entered by an operator through the I/O controller 
4, and the output unit 5 for outputting the pattern shape of the 
resist. 

With reference next to Figs. 10, 11 and 12, the method 
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for development profile simulation using the apparatus for 
development profile simulation shown in Fig. 9 according to 
the first embodiment of the present invention is described 
below. In Figs. 10 and 11, a method for optimizing development 
5 parameters "k" and "a" used for the method for development 
profile simulation is represented. 

In the case where the aerial image intensity is used as 
the optical intensities, the development parameters "k" and 
"a" are optimized through the processing in Fig. 5. However, in 

10 the case where the concentration of the reaction product is 
used as the optical intensities, the p hotoreaction parameters 
are also optimized by using a photoreaction model. Some 
measured photoreaction parameter is measurable, but other 
measured photoreaction parameter is not measurable. 

15 Therefore, not only development parameters but also 
photoreaction parameters are also optimized. The 
photoreaction by the exposure is different between a novolac 
resist and a chemically amplified resist. Therefore, the 
photoreaction model is different between the novolac resist 

20 and the chemically amplified resist. Further, -since the 
photoreaction of the chemically amplified resist is complicated, 
the photoreaction and the parameter affecting the pattern size 
is to be selected depending on the type of the resist. In Figs. 10 
and 11, acid generating efficiency, reaction rate constant, 

25 reaction order, and diffusion constant are selected as the 
photoreaction parameter. 
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(a) In step S201, N resist patterns are formed, and the N 
resist pattern sizes are measured. The measured pattern size 
Wi.o (i = 1, 2, 3, . . . , N) ofthe i*th resist pattern is stored in the 
data memory 2 shown in Fig. 9. The number of resist patterns 
(N) should be larger than the number of parameters to be 
optimized. When the number of the optimized parameters is x, 
N is set to x+1 or more. 

(b) In step S202, the measured initial dissolution rate 
R(co) of the developer used for pattern formation of the resist, 
and the measured dissolution rate R(c) during the development 
are measured. For example, a set of developers having a 
different alkaline concentration "c" are prepared, and the 
plural measured dissolution rate R(c) at the alkaline 
concentration "c" are measured. 

(c) In step S203, as described with reference to Figs. 1 to 
3, the ratio of the logarithm of the measured dissolution rate 
R(c) to the alkaline concentration c is calculated by using the 
measured dissolution rate R(c) as the measured value of the 
changing ratio of the dissolution rate. Thereafter, a table is 
formed containing the measured value of the changing ratio of 
the dissolution rate "g". The measured value of the changing 
ratio ofthe dissolution rate "g" is either completely a constant 
value, or differs in accordance with the alkaline concentration 
of the developer, the resist position in the thickness direction 
and the exposure dose. Further, in the case that it is possible 
to obtain the approximate function of the logarithm of the 
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measured dissolution rate R(c) relating to the alkaline 
concentration "c" of the developer, the measured value of the 
changing ratio of the dissolution rate "g" is obtained from the 
fitting function. Since the measured initial dissolution rate 
R(co) depends on the exposure dose, the measured value R(co) 
is obtained as the table or the fitting function. The obtained 
measured value of the changing ratio of the dissolution rate 
"g" is stored in the data memory 2. 

(d) In step S204, the exposure condition used for the 
actual measurement, the condition for the mask pattern, such 
as the pattern arrangement and the mask size, a threshold 
value 6 that is described below and used as a reference for 
optimizing the development parameter, and the number of 
parameter sets "M" to be employed are entered in The 
photoreaction product concentration calculator 31. The 
conditions may be entered by reading data from the data 
memory 2, or by an operator using the input unit 6. The 
"exposure condition" used for the actual measurement includes 
the condition such as NA and a and the exposure dose,. The 
"parameter set" is a combination of the development 
parameters "a", "k" and photoreaction parameter to be 
optimized. The parameter set containing the first parameter 
set to the M-th parameter set is prepared. 

(e) In step S205, The photoreaction product 
concentration calculator 31 calculates the aerial image 
intensity I(x, y) by using the conditions entered in step s204. 
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Thereafter, in step S206, The photoreaction product 
concentration calculator 31 reads the initial measured 
dissolution rate R(co) stored in the data memory and the ratio 
of the logarithm of the measured dissolution rate to the change 
of the alkaline concentration given by the equation (3). 

(f) In step S207, The photoreaction product 
concentration calculator 31 reads the n-th parameter set from 
the data memory 2. The n-th parameter set includes the 
photoreaction parameter, development parameter "oc n " and "k n ". 
The photoreaction parameter includes the acid generating 
efficiency C e ff,n, . . the reaction rate constant ki, n , k2,n, . . ., 
the reaction order p n , qn, . . ., and the diffusion constant Di,n, 
D2,n, . . For the first calculation, n is defined as 1. In step 
S208, The photoreaction product concentration calculator 31 
calculates the concentration M(x, y) of the photoreaction 
products using the aerial image intensity I(x, y) and the n*th 
parameter set. 

(g) In step S209, the spatial average value calculator 8 
calculates a spatial average value G(x, y) by the equation (6) 
with the calculated concentration M(x, y) of the photoreaction 
products. Here, I(x, y) in the equation (6) is replaced with M(x, 
y). In step S210, the dissolution rate calculator 10 calculates 
the calculated dissolution rate R(x, y) by using the spatial 
average value G(x, y) and the measured changing ratio of the 
dissolution rate of the resist. Thereafter, in step S211, the 
pattern shape calculator 11 estimates the resist pattern shape 
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by using the calculated dissolution rate R(x, y) and the 
concentration M(x, y) of the photoreaction products obtained in 
step S208. In step S212, the pattern shape calculator 11 
obtains a calculated pattern size Wi, n (i = 1, 2, 3, . . N) in a 
target position of the i~th resist pattern that corresponding to 
the measured pattern size Wi,o. 

(h) In step S213, the pattern shape calculator 11 judges 
the similarity between the calculated pattern size Wj >n and the 
measured pattern size Wi,o based on the threshold value 5. 
Specifically, if the sum value of the squared difference between 
the measured pattern size Wi >0 and the calculated pattern size 
Wi, n for the N resist patterns represented as Z| Wi, n - Wj,o I 2 is 
smaller than the squared threshold value 8 2 , step S214 is the 
next procedure. When Z| Wi, n * Wi,ol 2 is equal to or larger than 
the squared threshold value 8 2 , step S215 is the next procedure. 
It should be noted that 8 = (S | Wij - W i>0 I 2 ) 1/2 is employed as a 
threshold value for the first calculation (n=l). 

(i) In step S214, the pattern shape calculator 11 defines 
8 = (S | Wi, n " Wi,ol 2 ) 1/2 as a new threshold. In this case, the 
photoreaction parameter C e ff,n, ki, n , k2,n, Pn, q n> Di, n , D2,n and 
the development parameter "a" and "k" are stored in the data 
memory 2 as the temporary parameters. 

In step S215, the pattern shape calculator 11 determines 
whether "n" is equal to "M". When "n" is equal to "M" which is 
the number of the parameters set, the processing in Figs. 10 
and 11 is terminated. When "n" is not equal to "M", it is 
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assumed that a parameter set to be calculated remains and 
step S216 is the next procedure. 

(j) In step S216, the pattern shape calculator 11 
prepares the n+l-th parameter set, and performs the processes 
5 from steps S207 to S215 using the new parameter set. The 
processes from steps S207 to S215 are repeated until "n" is 
equal to "M'\ 

As described above, by searching for the parameter set 
minimizing the Z | Wi, n - W i)0 I 2 which is sum value of the 
10 squared difference between the measured pattern size Wi, 0 and 
the calculated pattern size Wi, n , the development parameters 
"a" and "k" and the photoreaction parameter are optimized. 

It should be noted that the photoreaction model and the 
photoreaction parameters are depend on the type of the resist 
15 and the appropriate reaction model is selected for each type of 
the resist. 

With reference next to Fig. 12, the method for 
development profile simulation using the optimized 
development parameters "a" and "k" and the photoreaction 
20 parameters according to the first embodiment is described 
below. 

(A) In step S221, the exposure condition for the resist 
pattern for which the size is to be predicted, the mask pattern 
conditions such as the pattern arrangement and the mask size, 
25 the measured initial dissolution rate R(eo), the measured 
changing ratio of the dissolution rate "g", and the optimized 
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development parameters "a" and "k", and the optimized 
photoreaction parameters in the data memory 2 are entered in 
the aerial image intensity calculator 7. The conditions may be 
entered by an operator using the input unit 6. The exposure 
condition includes the condition such as NA and a and the 
exposure dose. The development parameters "a" and "k" and 
the photoreaction parameters are the optimized parameters 
obtained through the processing in Figs. 10 and 11. 

(B) In step S222, the photoreaction product 
concentration calculator 31 calculates the aerial image 
intensity I(x, y) for a predetermined position (x, y) of the resist 
with a set condition such as exposure and the mask pattern. 
In step S223, The photoreaction product concentration 
calculator 31 calculates the photoreaction product 
concentration M(x,y) with the calculated aerial image 
intensity I(x, y) and the photoreaction parameters optimized 
through the processing in Figs. 10 and 11. 

(C) In step S224, the spatial average value calculator 8 
calculates the spatial average value G(x, y) of the aerial image 
intensity by using the equation (6) with the calculated 
concentration M(x, y) of the photoreaction products. Here, I(x, 
y) in the equation (6) is replaced with M(x, y). In step S224, 
the dissolution rate calculator 10 calculates the calculated 
dissolution rate R(x, y) for the predetermined position of the 
resist by using the average spatial value G(x, y), the measured 
initial dissolution rate R(co) , the measured changing ratio of 
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the dissolution rate "g", and the optimized development 
parameters "a" and "k". 

(D) In step S226, the pattern shape calculator 11 
predicts the pattern shape of the resist by using the calculated 
dissolution rate R(x, y) and the aerial image intensity I(x, y) 
obtained in step S223. Thereafter, in step S227, the pattern 
shape calculator 11 predicts the calculated pattern size W for 
the target position. 

In the above described method for development profile 
simulation, the dissolution rate in the predetermined position 
of the resist varies in accordance with the average spatial 
value of the optical intensities at the predetermined position 
of the resist. The measured changing ratio of the dissolution 
rate "g" used in the method for development profile simulation 
is defined in accordance with at least one of the exposure dose 
on the resist, the resist position in the thickness direction, and 
the alkaline concentration of the developer. Therefore, the 
changing ratio of the dissolution rate "g" a relating to the 
average spatial value of the photoreaction product 
concentration is also defined in accordance with at least one of 
the exposure dose on the resist, the resist position in the 
thickness direction, and the alkaline concentration of the 
developer. Consequently, the method for development profile 
simulation reflecting the actual behavior of the development is 
realized. Further, the method for development profile 
simulation according to the first embodiment makes it possible 
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to accurately predict the shape and the size of the resist 
pattern. 

In the method shown in Fig. 11, the threshold value 5 is 
updated in order to minimize Z | Wi, n * Wi,o I 2 as the result of 
the squared difference between the measured pattern size Wi,o 
and the calculated pattern size Wi, n for each pattern. However, 
so long that the S I Wi, n - Wi,o I 2 is minimized, a method is not 
limited to that shown in Fig. 11. 

Next, the results obtained by several types of simulation 
for various specific resist types are described below. Prior to 
describe the result by applying the method for development 
profile simulation in accordance with the first embodiment to a 
line pattern of the resist, the results obtained by a method 
disclosed in Japanese Patent Laid-Open Publication No. Hei 
10-64792 is described below. 

In the publication of Japanese Patent Laid-Open 
Publication No. Hei 10-64792, the changing ratio of the 
dissolution rate is proportional to the average spatial value of 
the intensity for a latent image or a photosensitizer 
concentration. Further, the changing ratio of the dissolution 
rate is assumed as a constant without actual measurement, 
and is optimized as one of the development parameters. Since 
the changing ratio of the dissolution rate is defined as being 
proportional to the average spatial value of the latent image 
intensity or the photosensitizer concentration, the equation 
(10) is used instead of the equation (7). 
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R(x, y) = R(co)[l + y{G(x, y)-G 0 >] (10) 

Under the condition, "y", "k" and "Go'' are defined as the 
5 development parameters and are optimized through the 
processing in Fig. 5. In this case, the changing ratio "g" is 
neither actually measured nor read from the data memory 2. 
A 1 = 1 line and space (L&S) pattern from 0.13 jam to 0.2 jam and 
the exposure condition, for example, the cyclic illumination 

10 with NA = 0.6, a = 0.75 and e = 2/3 are employed as the 
condition to optimize the parameters. Further, a resist of 300 
nm thicknesses on an anti -reflection film is used. 

However, when the "resist C" is used, the development 
parameters is not converged and consequently it is impossible 

15 to execute the optimizing method shown in Fig. 5. Therefore, 
the resist profile simulation is executed by using the 
development parameters "k" and "get" as the development 
parameters with the equation (7). Note that in this case, the 
changing ratio "g" is also neither actually measured nor read 

20 from the data memory 2. Thereafter, the parameters is 
converged and optimized. 

Next, the final pattern size is predicted in the case of 
changing space width of the line pattern of 0.13 jam using the 
resist C. The result shows that the line width is formed 

25 smaller at the position where the ratio of the line width to the 
space width is 1-2 or 1'3 than at the position where ratio of the 
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line width to the space width is 1-1 when the line width of the 
mask is constant and the space width is varied. Furthermore, 
the line width increases as the space width becomes larger. 
Such a behavior of development is common for almost all types 
of the resist. When the simulation is executed by using the 
development parameter optimized by actually using the resist 
C to obtain the final line size, not only such a behavior of 
development is reproduced, but also a predicted pattern size is 
similar to the measured pattern size. 

However, actual measurement of the changing ratio of 

the dissolution rate "g" of the resist C represented by the 
equation (3) reveals that the resist C has character the 
changing ratio of the dissolution rate "g" does not depend on 
the concentration of the developer, the exposure dose and the 
resist position in the thickness direction. The measured 
changing ratio of the dissolution rate "g" is also substantially 
equal to the optimal value obtained by the simulation. 
Therefore, using the equations (3) and (7) may be appropriate. 
Further, in case of the resist having the substantially constant 
changing ratio of the dissolution rate "g" relating to the 
concentration of the developer, the exposure dose and the 
resist position in the thickness direction, it is possible to 
optimize the changing ratio "g" as the one of the development 
parameter obtained through the processing in Fig. 5 without 
actual measurement of the changing ratio "g". 

However, the above method is not available for all 
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conditions. According to the method, it is possible to optimize 
the development parameters "a" and "k" for the "resist A" 
shown in Fig. 2. The resist A has the changing ratio of the 
dissolution rate "g" varying at the resist position in the 
5 thickness direction. However, it is impossible to accurately 
predict the final pattern size with the optimized parameter 
when the space width with the line width of 0.13 \xm is changed. 
Consequently, though the above described method reproduces 
the behavior of the line width when the space width is changed, 

10 the predicted pattern size is not equal to the measured pattern 
size. Therefore, when the value of the changing ratio "g" is 
defined for each of the upper portion, the interior and the 
lower portion of the resist as described above, the predicted 
pattern size reproduces the measured pattern size. 

15 In addition, it is impossible to converge the development 

parameters by the above method when the "resist B" shown in 
Figs. 3A and 3B having the changing ratio of the dissolution 
rate "g" drastically changing in accordance with the alkaline 
concentration of the developer, the exposure dose and the 

20 resist position in the thickness direction. It is assumed that 
the development parameters is not converged because the 
changing ratio of the dissolution rate "g" of the resist B 
drastically changes as shown in Figs. 3A and 3B when the ratio 
is measured. Therefore, after the changing ratio of the 

25 dissolution rate "g" is measured, a table is formed containing 
the measured changing ratio "g" for each conditions such as 
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the concentration of the developer, the exposure dose and the 
resist position in the thickness direction. Thereafter, the 
development parameters "a" and "k" are optimized by using the 
value in the table of the changing ratio of the dissolution rate 
5 "g". As the result, the convergence value for the development 
parameters is obtained. Thereafter, the final pattern size is 
predicted by changing the space width with the line pattern of 
0.13 jam. Consequently, the behavior of the line width when 
the space width changes is reproduced. Further, the 

10 predicted pattern size is substantially equal to the measured 
pattern size. Therefore, forming the table containing the 
actually measured changing ratio of the dissolution rate "g" 
and executing the optimization of the parameters and the 
pattern resist profile simulation provide precise prediction of 

15 the pattern size. 

In the first embodiment, the ratio of the logarithm of the 
measured dissolution rate R(c) to the alkaline concentration is 
employed as the changing ratio of the dissolution rate of the 
resist as shown by the equation (3). However, the equation 

20 representing the changing ratio of the dissolution rate is not 
limited to the equation (3). For example, as shown by the 
equation (2), when the logarithm of the measured dissolution 
rate R(c) is proportional to the logarithm of the alkaline 
concentration, the changing ratio of the dissolution rate "g" 

25 may be given by the equation (ll). 
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a{log R(c)} /d{log c} = g (11) 



Further, the dissolved amount of the resist around the 
resist pattern is obtained as the average spatial value G(x, y) 
of the aerial image intensity I(x, y) in the first embodiment as 
shown by the equation (6). However, the average spatial 
value G(x, y) is not limited to the description. When the 
positive resist is used, an approximation function H(x, y) given 
by the equation (12) may also be available. 

H(x,y) 

=JJdx , dy , I(x,y)f(a)exp[{-(x-x') 2 +(y-y , ) 2 }/k2] (12) 

Here, f(a) is 1 when I(x, y) > a, while f(a) is 0 when I(x, y) 
< a. This is based on the following reason. Since the exposure 
dose sufficient to form the resist pattern makes the dissolution 
rate of the resist very high such as several thousands nm/sec, 
an area exposed in a certain amount or more is dissolved in the 
initial stage of development. Thereafter, the development of an 
area for determining the size of the pattern edge portion is 
initiated. Therefore, it is assumed that, as the high exposure 
area of the resist dissolves, the alkaline concentration is 
changed, and the developer having a reduced alkaline 
concentration makes the pattern edge portion dissolve. 
Consequently, so long as the area where the aerial image 
intensity I(x, y) is equal to or larger than a constant value (a), 
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the other area where the aerial image intensity I(x, y) is less 
than the constant value (a) can be neglected. As the result, it 
is possible to reduce the calculation time. On the other hand, 
in case of the negative resist, f(a) is 0 when I(x, y) > a, while 
f(a) is 1 when I(x, y) < a in the equation (12). 

When the approximation function H(x, y) given by the 
equation (12) is used, the concentration M(x, y) of the 
photore action products reacted by the exposure and the PEB is 
also available instead of I(x, y). In this case, it is possible to 
optimize the constant value (a) in the equation (12) for the 
simulation as a development parameter as well as a and k. 

As described above, the method for development profile 
simulation in accordance with the first embodiment makes it 
possible to make a phenomenon model such that the large area 
dissolves at the initial development stage, the alkaline 
concentration of the developer is changed, and the developer 
makes the low exposure area dissolve, so that the final size is 
determined. Therefore, method for development profile 
simulation according to the first embodiment makes it possible 
to reduce the calculation time without losing the accuracy. 

Especially, the spatial average value calculated in the 
area where the aerial intensities or the concentration of the 
photoreaction products is either equal to or larger than the 
constant value or less than the constant value is used. 
Therefore, it is possible to provide the phenomenon model that 
the large area dissolving at the initial development stage 
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changes the alkaline concentration of the developer and the 
developer makes the low exposure area dissolve and the final 
pattern size is determined. Consequently, it is possible to 
reduce the calculation time while the accuracy is maintained. 

The method for development profile simulation 
according to the first embodiment of the present invention is 
capable of being expressed as descriptions of a series of 
processing or instructions for a computer system. Therefore, 
the method for development profile simulation is capable of 
being formed as a computer program product to execute 
multiple functions of the CPU in the computer system. "The 
computer program product" includes, for example, various 
writable mediums and storage devices incorporated or 
connected to the computer system. The storage device is 
capable of being expressed as the program storage unit 3 
shown in Fig. 4. Writing the computer program into the 
program storage unit 3 is also available. Consequently, the 
computer program of the method for development profile 
simulation is executed by the processing unit 1 in accordance 
with predetermined procedures. The writable mediums include 
a memory device, a magnetic disc, an optical disc and any 

devices that record computer programs. 

The computer system may have a flexible disk drive and 

a CD-ROM drive in front. The computer program product as the 

flexible disk or the CD-ROM is inserted into each disk drive 

and executes the method for development profile simulation. A 
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ROM as a semiconductor memory for example used for game 
cassette and a cassette tape as a magnetic tape are also 
available for the computer program product by connecting a 
drive unit to the computer system. 

(Comparison with An Actual Value) 

As an example that the method for development profile 
simulation according to the first embodiment is applied, the 
method for predicting pattern size for a chip having the same 
pattern arranged in array is described below. In Fig. 7, a 
layout pattern includes a non-exposure resist portion 26 and a 
hole pattern 25 surrounded by the non-exposure resist portion 
26. The hole pattern 25 has a plural hole patterns arranged 
in array. The positive resist of 300 nm thickness coated on 
the anti-reflection film on the wafer. The diameter of the hole 
in the hole pattern is 0.2 ^m and the ratio of the diameter of 
the hole to the space is 1-1. While the developer is supplied 
through a linear nozzle, the wafer is rotated to apply the 
developer uniformly. The annular illumination with NA = 
0.68, a = 0.75 and s = 2/3 are used as the exposure condition. 

With reference next to Figs. 8A and 8B, "Pa" and "Pb" 
along the horizontal axes correspond to the position of the hole 
pattern shown in Fig. 7. As shown in Fig. 8A, for the 
measured resist pattern, the hole diameter in peripheral 
position of the hole pattern 25 such as Pa and Pb is larger than 
the hole diameter in central region where the hole patterns are 
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concentrated. However, when the resist profile simulation is 
executed by only using the aerial image intensity as shown in 
Fig. 8B, the hole diameter in Pa and Pb is smaller than the 
hole diameter in the central region where the patterns are 
concentrated. 

The large area adjacent to the points Pa and Pb is 
non-exposure area. When the developer is supplied to the 
rotated wafer through the linear nozzle, there is almost no 
movement of the developer and only the diffusion of the 
dissolved material occurs. Therefore, since the consumption 
of alkaline in the developer around Pa and Pb is smaller than 
the consumption of alkaline around the central region of the 
hole pattern portion 25, the regions around Pa and Pb are 
developed with the developer of high alkaline concentration. 
Therefore, it is assumed that the final hole diameter is 
increased. 

Therefore, the changing ratio of the dissolution rate "g" 
of the resist D is measured and the method for development 
profile simulation according to the first embodiment is 
executed. The development parameters are optimized by 
changing the ratio of the hole diameter to the space while 
maintaining the hole diameter of 0.2 ^m. As the result, it is 
possible to reproduce the measured hole diameter shown in Fig 
8A by the simulation. Consequently, it reveals that the 
method for development profile simulation according to the 
first embodiment predicts the pattern size difference in the 
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chip even when the density of the chip patterns is distributed. 

(Second Embodiment) 

A method for predicting the size difference due to flow of 
5 developer with profile simulation is described in a second 
embodiment. When the developer is being supplied, or after 
the developer puddle is formed, flow of the developer may 
occur depending on the process condition. In this case, the 
flow velocity of the developer is higher than the diffusion rate 

10 of the dissolved reaction product or the alkaline. With the 
normal exposure dose, the dissolution rate on the large area is 
very high such as several thousands nm/sec. Further, the area 
exposed certain exposure dose or more dissolves in a moment. 
Therefore, when the flow of the developer occurs, a target 

15 pattern is affected by the alkaline concentration of the 
developer flowing into the pattern and the dissolution rate is 
changed. 

With reference now to Fig. 13, the developer 23 is 
supplied to the resist 24 coated on the substrate uniformly. 

20 The developer 23 flows at a constant velocity from reference 
position A(ct,p) to a target position X(x, y). At time "t", the 
developer 23 flows from the reference position A to the target 
position X where a pattern shape is calculated. At the 
reference position A, since the developer 23 reacts with the 

25 resist 24, the alkaline concentration co of the supplied 
developer 23 is reduced to the alkaline concentration c. 
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Therefore, the alkaline concentration c of the developer 23 at 
the target position X(x, y) depends on the dissolved resist 
amount around the reference position A(a,p). Thus, the 
reference position A(a,p) on the resist 24 for calculating the 
aerial image intensity and the average spatial value of the 
intensity is defined in different position to the target position 
X(x, y) on the resist 24 for predicting the pattern shape. 
Further, the reference position A(a,p) for calculating the 
aerial image intensity and the average spatial value of the 
intensity be moved in accordance with the flow velocity of the 
developer 23 and the development time. Thereafter, reference 
position optical intensities I(ct,p) and the average spatial 
value G of the intensity are calculated at the reference 
position A(a,p). The reference position A(a,p) is calculated as 
a function (a (t), p (t)) for time t. 

For a positive resist, the area exposed certain amount or 
more dissolves in a moment. The dissolved resist amount 
around the reference position A during a constant time of the 
initial development stage is calculated under assumption that 
the dissolved resist amount is proportional tt) the average 
spatial value of the aerial image intensity. H(x, y) given by 
the equation (12) is used as a concentration modulation 
function instead of G(ct,p) given by the equation (6). 
Therefore, the equation (14) where variables x and y in the 
equation (12) are replaced with a and p is used. 
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H(a,p) 

= JJdx f dy'I(a,p)f(a)exp[{-(a-x , )2 + (p-y , ) 2 }/k 2 ] (14) 



Here, f(a) is 1 when I(oc,p)> a, while f(a) is 0 when 
5 I(a,p) < a. 

On the other hand, in the case of the negative resist, 
since the area exposed less~than the certain amount is 
dissolved, f(a) is 0 when I(a,p) > a, while f(a) is 1 when I(ot,p) < 
a in the equation (14). 

10 Similar to the first embodiment, the ratio of the 

logarithm of the measured dissolution rate R(c) of the resist to 
the alkaline concentration "c" of the developer is given by the 
equation (3) where the logarithm of the measured dissolution 
rate R(c) of the resist to the alkaline concentration "c" of the 

15 developer. When co is given as the initial alkaline 
concentration of the supplied developer, the measured 
dissolution rate R(c) is obtained by the equation (4). 
Therefore, the calculated dissolution rate R(x, y) is obtained 
by the equation (15). For each time "t", the calculated 

20 dissolution rate R(x, y) is given by the equation (15). 

log R(x,y) = log R(co) + gTiH(a (t), p (t)) (15) 

The development parameters are optimized through the 
25 processing in Fig. 5. In this case, the development parameter 
includes a constant exposure dose (a) and a development 
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parameter r\ instead of the development parameter a. The 
measured pattern size used for the optimization may be actual 
measurement under the condition where the development is 
performed with flowing developer. The actual measurement 
under the condition where the flow of the developer does not 
occur is also used as an alternative. 

Subsequently, the profile simulation is executed to 
predict the pattern size by using the optimized parameters and 
the method shown in Fig. 6. In this case, the aerial image 
intensity I(x, y) is replaced with the reference position optical 
intensities I(a,p). H(a,p) given by the equation (14) is 
employed as the average spatial value. The calculated 
dissolution rate R(x, y) is obtained by the equation (15). 
Further, the aerial image intensity in the target position X and 
the calculated dissolution rate obtained by the equation (15) 
are employed to predict the pattern shape. 

The predicted pattern size obtained by the profile 
simulation and the measured pattern size are compared by 
using the mask pattern same as the hole pattern shown in Fig. 
7. With reference to Fig. 14, the mask pattern according to 
the second embodiment includes an non-exposure resist area 
26 and a hole pattern 25 similar to the hole pattern shown in 
Fig. 7. In the hole pattern 25, a plural hole patterns are 
regularly arranged in array. Further, a scan development 
method which supplies the developer by scanning a linear 
nozzle from one end to the other end of the wafer in a nozzle 
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scanning direction is shown in Fig. 14 is used. Since the scan 
development method does not rotate the wafer, the scan 
development method is an effective development method for a 
wafer having a large diameter, a large photo mask and a liquid 
crystal display. However, while the flow of the developer is 
substantially negligible in the case where the wafer is rotated, 
the supplied developer gradually moves and the uniformity is 
deteriorated in the scan development method. The direction 
of the flow of the developer is opposite at 180 degrees to the 
nozzle scanning direction in Fig. 14. The linear nozzle is 
moved from the left to the right in Fig. 14. In this case, after 
the puddle is formed, the supplied developer moves from the 
right to the left in Fig. 14. Therefore, the developer 
maintaining the alkaline concentration reaches the position Pb 
and the hole diameter becomes large. 

With reference to Fig. 15A, the hole diameter in the 
position Pa is slightly larger than the diameter around the 
center of the hole pattern portion 25. The phenomenon is due 
to the affect of existing of the adjacent non-dissolved area at 
the initial development stage. Meanwhile, since the flow of 
the developer is not considered in the method for development 
profile simulation according to the first embodiment, th^ 
distribution of the diameter is obtained as shown in Fig. 8A 
and the measured diameter shown in Fig. 15A that needs to 
consider the flow of the developer is not predicted. Therefore, 
the profile simulation is executed by using H(a,p) given by the 
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equation (14) and changing the reference position along with 
the time to calculate the concentration. As the result, the 
hole size distribution is substantially equal to actual 
measurement shown in Fig. 15A and the calculated pattern 
5 size is also substantially equal to the measured pattern size. 
As described above, when the flow of the developer occurs, it is 
also possible to predict the change in the pattern size due to 
the flow of the developer. 

Consequently, in the case where the developer flows 
10 during the development process, it is possible to predict the 
varied pattern size in the target position under condition in 
which the alkaline concentration of the developer is reduced 
because a photosensitive resin dissolves in the reference 
position. 

15 In the second embodiment, H(a,p) is used as the 

concentration modulation function. However, the equation 
(6) may be used by replacing the variables x and y with a and P 
Therefore, G(ct,p) given by the equation (16) is used as an 
alternative . 

20 

G(ot, p) 

= JJdx'dy' I(cc,p) exp [{-(a-x') 2 +{p -yO^/k 2 ] (16) 

25 (Third Embodiment) 

It is possible to apply the method for development 
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profile simulation according to the first and second 
embodiment to optimization of mask pattern data in order to 
reduce projected pattern size variation. A method for mask 
pattern data correction using a predicted resist pattern shape 
in accordance with a third embodiment is described below. 

With reference now to Fig. 16, an apparatus for mask 
pattern data correction according to the third embodiment of 
the present invention includes a processing unit 1 having a 
function for performing the method for mask pattern data 
correction, a data memory 2 and a program memory 3 
connected to the processing unit 1. The data memory 2 stores 
various data, such as an exposure condition, mask pattern data 
and a development parameter 12 which are used for mask 
pattern data correction. The program memory 3 stores 
various types of application software, middleware and an 
operating system that are used for mask pattern data 
correction. 

The processing unit 1 includes an aerial image intensity 
calculator 7 configured to calculate an aerial image intensity 
of a resist, a spatial average value calculator 8 configured to 
calculate a spatial average value of the aerial image intensity, 
a changing ratio of the dissolution rate reader 9 configured to 
read the measured changing ratio of the dissolution rate of the 
resist relating to an alkaline concentration which depends on 
at least one of the exposure dose, the resist position in the 
thickness direction, and the alkaline concentration of the 
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resist, a pattern reader 41 configured to read mask pattern 
data and designed pattern data in the resist, a dissolution rate 
calculator 10 configured to obtain a calculated dissolution rate 
using the spatial average value and the measured changing 
ratio of the dissolution rate, and a pattern shape calculator 11 
configured to predict the pattern shape of the resist by using 
the calculated dissolution rate. The processing unit 1 is 
connected to an input unit 6 for receiving data and 
instructions entered by an operator through an I/O controller 4, 
and an output unit 5 for outputting the pattern shape of the 
resist. 

With reference next to Fig. 17, the method for mask 
pattern data correction using the apparatus for mask pattern 
data correction is described below. 

(a) In step S321, the designed pattern size Wi,o (i = 1, 2, 
3, . . ., N) in target position of the designed pattern data of the 
resist, size allowance 5, exposure condition, mask pattern 
condition including pattern arrangement and mask size 
included in the mask pattern data are entered in the 
processing unit 1. The development parameters "a" and "k", 
measured dissolution rate R(co), and the measured value of the 
changing ratio of the dissolution rate "g" are also entered in 
the processing unit 1. Here, a plural mask pattern set having a 
different pattern arrangement and size are prepared. 

(b) In step S322, the mask pattern size Mi.o (i = 1, 2, 
3, . . ., N) is entered in the aerial image intensity calculator 7. 



Thereafter, in step S323, the aerial image intensity calculator 
7 calculates the aerial image intensity I(x, y) with the mask 
pattern size Mi.o as a variable. 

(c) In step 324, the spatial average value calculator 8 
calculates a spatial average value G(x, y) of the aerial image 
intensity by the equation (6) with the calculated aerial image 
intensity I(x, y). In step S325, the dissolution rate calculator 
10 calculates a calculated dissolution rate R(x, y) by using the 
spatial average value G(x, y), the measured initial dissolution 
rate R(co), the measured changing ratio of the dissolution rate 
"g", and the development parameter "a" and "k" in the 
reference position of the resist. Thereafter, in step S326, the 
pattern shape calculator 11 estimates the resist pattern shape 
by using the calculated dissolution rate R(x, y) and the aerial 
image intensity I(x, y) obtained in step S323. In step S327, the 
pattern shape calculator 11 calculates the calculated pattern 
size Wi, n in the target position. 

(d) In step S328, the pattern shape calculator 11 judges 
the similarity between the calculated pattern size Wi, n and the 
designed pattern size Wi,o based on the size allowance 6. 
Specifically, if the sum value of the squared difference between 
t^ie designed pattern size Wi,o and the calculated pattern size 
Wi, n for 1 st to N-th target position represented as Z| Wi >n " 
Wi,o I 2 is smaller than the squared value of the size allowance 
5 2 , the procedure shown in Fig. 17 is terminated. When I I Wi, n 
" Wi,o | 2 is equal to or larger than the squared value of the size 
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allowance 8 2 , it is determined that the mask pattern size Mi.o 
entered in step S322 does not reproduce the designed pattern 
size of the designed resist pattern data on the resist. In this 
case, step S329 is the next procedure. 

(e) In step S329, the pattern shape calculator 11 changes 
mask pattern set. Specifically, the mask pattern size Mi,„ is 
replaced with Mi, n + i. Thereafter, the processes from steps S322 
to S328 are repeated with the replaced mask pattern size. In 
step S330, the mask pattern data is corrected by using the 
optimum mask pattern size Mi, n . 

As described above, by searching for the parameter set of 
mask pattern minimizing the I I Wi >n - Wi,o I 2 which is sum 
value of the squared difference between the designed pattern 
size Wj,o and the calculated pattern size Wi, n , it is possible to 
optimize the mask pattern size and fabricate the optimized 
mask pattern. It should be noted that the concentration M(x, y) 
of photoreaction products is an alternative for the aerial image 
intensity I(x, y) in the method for mask pattern data 
correction. 

Since the method for mask pattern data correction 
makes it possible to reduce the pattern size variation 
depending on the development process, the desirable resist 
pattern shape is accurately formed on the resist. 

With reference again to Fig. 7, it is possible to reduce 
the difference between the hole diameter in the central region 
and one in Pa and Pb with the method for mask pattern data 
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correction. In this case, the hole diameter in the central region 
is calculated by the profile simulation and the calculated hole 
diameter is used as the designed pattern size Wi,o. After the 
calculated pattern size Wi, n of each hole pattern located on 
peripheral region is obtained, the mask pattern size Mi, n 
corresponding to the hole pattern located on peripheral region 
is optimized. By using the optimized mask pattern data, the 
optimized mask is fabricated. Without optimizing the mask 
pattern, the hole diameter in the peripheral region is larger 
than the hole diameter in the central area, as shown in Fig. 
15A. However, by using the method for mask pattern data 
correction according to the third embodiment, it is possible to 
obtain uniform hole diameter independent of the region, as 
shown in Fig. 15B. 

In the method shown in Fig. 17, the mask pattern size 
Mi, n is optimized until I | Wi, n - W i>0 I 2 is smaller than the 
squared value of size allowance 5 2 . However, so long that the 
E|Wi, n • Wi,o I 2 is minimized, a method is not limited to that 
shown in Fig. 17. As the method shown in Figs. 5 and 11, 
when I|Wi, n * Wi,o I 2 is smaller than the squared size 
allowance value 5 2 , defining 5 = (Z | Wi, n - Wi.ol 2 ) 172 as a new 
size allowance and repeating the judgment of the similarity 
between the calculated pattern size Wi, n and the designed 
pattern size Wi,o with all prepared mask pattern Mi, n is an 
alternative . 

In the first embodiment, the method for development 
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profile simulation to predict the difference between the hole 
diameter in the central region and one in the peripheral region 
Pa, Pb shown in Fig. 7 and Fig. 8A is described. The method for 
mask pattern data correction according to the third 
embodiment makes it possible to reduce the difference by 
changing the mask pattern size. 

(OTHER EMBODIMENTS) 

Although the invention has been described above by 
reference to the embodiment of the present invention, the 
present invention is not limited to the embodiment described 
above. Modifications and variations of the embodiment 
described above will occur to those skilled in the art, in the 
light of the above teachings. 

The method for development profile simulation 
according to the first and second embodiment makes it possible 
to predict the pattern size variation occurred in the 
development process. Further, since the method for 
development profile simulation predicts the pattern size in the 
target position, it is possible to predict the distribution of the 
variation of the pattern size on a wafer or an LSI chip. 
Therefore, it is possible to determine whether the variation of 
the predicted pattern size is allowable or not. 

Further, the method for development profile simulation 
according to the first and second embodiment is useful for 
selecting the condition of lithography process prior to actual 
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fabrication. In this case, the parameter set defined in step 
S104 shown in Fig. 5 includes not only the mask pattern size 
but also the exposure condition such as NA, a and the 
illumination condition. Using the parameter of the exposure 
condition increases the precision of the predicted pattern size. 

The method for development profile simulation 
according to the first and second embodiment is also useful for 
selecting the resist. The method for development profile 
simulation makes it possible to select the optimum resist from 
the plural types of resist depending on the type of the 
semiconductor device or the type of a lamination layer. By 
preparing the table containing the measured changing ratio of 
the dissolution rate "g" for the plural type of the resist, it is 
possible to select the optimum resist. Though actually 
manufacturing masks over and over requires a great deal of 
labor, time and costs, the method for development profile 
simulation selects the optimum resist worth using for 
manufacturing. Therefore, it is possible to reduce development 
time. 

It should be noted that the term of "optical image" means 
the projection image formed on wafer by light in narrow sense. 
However, the optical image means both the projection image 
and a latent image in the broad sense. In the first and second 
embodiment, the term of the optical image means the 
projection image. However, the latent image is available 
instead of the projection image. Especially, the term of the 
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optical image in the following claims includes both the 
projection image and the latent image. 

As described above, the present invention includes many 
variations of embodiments. Therefore, the scope of the 
5 invention is defined with reference to the following claims. 
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